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. Abstract 

, Within the Chiral Quark Model we have calculated the effect of non-factorizable (G^) con- 

' tribution to the Bk parameter. When all diagrams are summed, a vanishing result is obtained. 



> ; 1 Introduction 

The low energy effects of X° — K mixing has been a subject of a lot of discussion in the literature. 
See for example and references therein. Its amplitude is parametrized by the Bk parameter 

, and is conventionally defined by 

A 



oo 



(X"|0(A5 = 2)\K") = -fK^j^Bj, (1) 

The vacuum saturation approach accomodates only a factorizable contribution which yields Bk = 
pi i' 3/4 in the large Nc limit. By adding the factorizable next to leading order contribution 1/Nc increases 

^ • Bk to 1. Lattice calculations seem to indicate that the physical Bk is numerically around 0.7. 

Non-perturbative and non-factorizable effects can be taken into account by dressing the Feynman 
diagrams with soft external gluons in all possible ways. In doing this one makes for instance the 
■ identification GJ^j.G"'^'' (G^) where G^i, is the gluon tensor of dimension two. The coefficient of 

the (G^) term has been calculated by Pich and de Rafael They got a significant contribution to 
Bk- Numerically, their result change very much the next to leading l/Nc contribution. However, 
due to large uncertainties in the experimental value of the condensate (G^), Bk is determined to lie 
between 0.3 and 0.6. The purpose of this paper is not to find a better numerical determination of 
Bk, but rather to find a more general expression of Bk in terms of non-perturbative (G"^) effects. 

In this paper we use the Chiral Quark Model with a constituent quark mass M to calculate the 
amplitude. This model introduces interaction terms between quarks and mesons. This means that 
hadronic matrix elements can be calculated as loop diagrams. 

When calculating fermion propagators in an external field we use the point gauge prescription 
in order to get results which directly contain gauge invariant quantities. 

There are two independent condensates with dimension six, namely (G'^) and (j^). Since we are 
working in the heavy quark limit, (j^) will be given as an expansion in gluon condensates of higher 
order and inverse proportional to the heavy quark mass M ^ i.e. 

To consider {j'^) terms we have to consider direct (C^) terms as well. Therefore, our object of 
investigation is the (G^) terms only. 
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2 The Chiral Quark Model 



In the Chiral Quark Model, chiral symmetry breaking is taken into account by adding an extra term 
to the standard QCD Lagrangian: 

where 

llQCD^q{ll-D-M^)q-\G'^^>^''G%. (4) 

and the L-^ term is only taken into account for energy/momenta of order A^^ and below. The q- fields 
denotes only u, d, s quarks. The chiral symmetry breaking scale is = 0.83 GeV and 

L^^-M{q^Y.qn + q^Y.\i:). (5) 

The prefactor M is interpreted as the constituent quark mass and is of the order 200 — 300 MeV . 
The pseudoscalar Goldstone-octet fields tt" are contained in 

E = exp(z^AV//) (6) 

a 

In the limit where M.q 0, the model is SU{3)l ® SU{3)r invariant. 

The chiral symmetry is thus realized through the existence of pseudoscalar Goldstone bosons. 

This model connects mesons and quarks and one can calculate hadronic matrix elements as loop 
diagrams. By expanding S one finds the coupling between a quark density and a meson. For the 
qqmeson vertex, the relevant coupling constant is oc M"f^/ f where / is associated with the physical 
pion decay constant and the proportionality factor is given by SU (3) symmetry. In the following, 
mesons are drawn by dotted lines while quarks are drawn with solid lines. 

In the SU{3) limit without gluonic condensates the pion decay constant is defined by 

f - (7) 
~ 47r2/ 

« " -O^ "(vertex : 7''L) (8) 

With a cut-off A, the explicit expression for = ln{A'^ / M'^) + ■ ■ ■ while in dimensional regularization 
/tt = {2 — D/2)~^ + ■ ■ ■. A is of the order of the chiral symmetry breaking scale Aj^ . Although/^ and 
/ appear difi'erently in the Chiral Quark Model, they are for practical purposes numerically equal 
(/^ = 93.3 MeV). 

3 Diagrammatic Calculation of — 

The lowest order local operator, of dimension six, contributing to — is 

©(AS* = 2) = (dj^Ls)(d-f''Ls). (9) 
This operator gives rise to two diagrams pictorially given as 

(XO|0(A5 = 2)|if) = " -OxO- + -<X>- " (10) 

The first diagram is cx while the second diagram is colour suppressed and is oc iVc. 

There are also other contributions than eq.^ with derivatives acting on the quantum fields. 
They are generated from the Siamese penguin operator |^], Q or similar non-local operators. 
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In some calculations it is customary to Fierz transform the operator O. We can take some 
advantages of the transformed operator as some of the generated diagrams cancel. Consider the 
Fierz transformed operator: 

Of - ■^(d^^Ls){d^^Ls) + 2{d^^Les){d^''Les) (11) 

The matrix element can now be written symbolically: 

{K"\Op\lf) = i-{"-OxO- + -<X>- {vertex -.^L) 

+2|" -OxO- + -OI>- '■'] {vertex -.-/^W) (12) 

The Fierz transformed operator gives rise to four diagrams. The symbolic descriptions of O and 
Of are of course identical. The first term in eq.|l^ corresponds to the second and the fourth term 
of eq.|lj. Similarly, the second term of eq.|l^ corresponds to the first and the third term of eqjl^. 
Examples below show that this is indeed the case whether effects of gluon condensates are included 
or not. One can in fact choose to consider only " factorizable" contributions with different vertices 
and prefactors (first and third term of eq.|l^ and first term of eq.p^. 

If the external field is turned off, one sees from the colour structure that the first term in eq.p^ 
corresponds to the second and fourth term in eq.|l2|: 

= -^Nc + 2rr(rr). (13) 

Similarly, the second term in eq.^ corresponds to the first term (and a vanishing third term) in 
eq. 
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In order to calculate non-perturbative contributions to Bk, we need fermion propagators in an 
external gluon field. A pedagogical presentation for this calculation is presented by Novikov et al. 
in j^. In the point gauge the gluon field can be expanded in terms of the gluon tensor G^q: 

A^{x) = ^x^Gp^iO) + ^x^x^[D^Gp^m + ^x^x''x^'[D^D,Gp^m + ■■■. (14) 

When calculating the fermion propagator in an external gluon field we see that we get contributions 
from "G", "DG" and "DDG" terms. In order to obtain gluon condensate contributions one has 
interpret products of gluon tensors as vacuum averages. An averaging over gluon indices yields for 
example the substitution 

gab 

Ga/jG^i/ ^ -^(G^) (50^5/31^ - QauOflu)- (15) 

This is how non-perturbative effects are interpreted from a perturbative prescription. Similar expres- 
sions as in eq.^, albeit more complicated, can be obtained with three gluons and terms including 
derivatives. They can be derived by making use of the equations of motion 

D^Gl, - -gsj:. (16) 

the commutation relation 

[D^,,D,]^~igsG^,, (17) 

the Bianchi identity and translational invariance. This will enable us to separate the (G'^) contribu- 
tion from (j^). We find the relations 

{Q\{D^G,,r{D'^GP^r\Q) = 2gl{f)-2g,{G') (18) 
{Q\{D^Gp^Y{DPG^^^r\Q) = g^Af)-9s{G') (19) 
{0\{D^G'^n''iDpGP,m = g'sif) (20) 
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See also reference |^ . 

When calculating the off-diagonal — K amplitude with external gluons, one can take ad- 
vantage of the diagrams already included in the expression for A more refined version of eq.0 
including gluon condensates can be written 



f 

J TT 



NeM"^ 

47r2 



,0C. 



— G') 



1 



(21) 



where the dots denotes higher order condensates. During the calculation we need some useful trace 
expressions of the SU (3) generators: 

(22) 



Tr{et^) = -(5''^ 



1 



and 



Trim") = -{d""" + ir'"') 



ANr 8 



(23) 



(24) 



When dressing the "blobs" from the Op operator with two external gluons ending in vacuum, 
one can realize that the two diagrams having only one single colour line, i.e. the second and fourth 
term of eq.|l^, indeed cancels in the two gluon case (one gluon on each loop) as seen from the colour 
structure 

(25) 



In addition, the first term in eq.^ vanishes due to traceless colour matrices. Thus, we only have to 
consider the diagrams coming from the third term in eq.|l^ with external gluons attached. This is a 
non-factorizable contribution due to the vertex structure 7^Lt°. Non-perturbative gluons to order 
(^G^) contributes to Bk as 
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IT 



2\1 



32^2/4 



(26) 



This is in accordance with The terms with two gluons on each loop are only contributing to /^r 
and therefore bring no effect into Bk- One should also have in mind that eq.^is further modified 
by meson loops . 

In the case of three gluons, one again realizes that the second and fourth term in eq.|l^ cancels 
in the case of two external gluons on one loop and one gluon on the other. The cancellation is again 
seen by looking at the colour structure. One encounters the sum 



—Tr{eth') + 2Tr{et^tW) = 



(27) 



where we have used the commutator and completeness relation of two SU (3) generators and applied 
the trace expressions for up to four SU{?>) generators. The diagrams with all three gluons attached 
to one loop are already included in the physical pion decay constant. Thus, we are left only with 
the third term in eq.|l| which we dress with gluons. Due to the different "G", "DG" and "DDG" 
parts in eq.^ there are totally 20 diagrams to calculate. 

A useful tool for the calculation is to use a software package for algebraic manipulation. We 
use the Form program, invented by J. Vermaseren [lll| , which is well suited for our purpose. By 
summing up all diagrams, we find zero contribution of the (G'^) terms. This result is a generalization 
of a result obtained in the paper of W. Hubschmid and S. Mallik [|l^ where they show that there are 
no contribution of (G'^) to two-point fimctions. A more complicated structure is encountered in our 
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case when calculating the parameter and their result is therefore not automatically applicable 
in our case. 

As a cross-check of our result, we also performed the calculation with just the original operator 
O. This was a more time consuming operation, but the same vanishing result was obtained. 

We have also calculated the contributions proportional to (j^) and obtained a non-vanishing 
constant term (constant with respect to the kaon momentum) which is not acceptable. However, 
(C^) type contributions should be included as well according to eq.|^. To complete this part of the 
calculation would require a major effort. 

4 Conclusion 

We have calculated the non-factorizable < G"^ > contribution to the Bk parameter. A vanishing 
result is obtained as did for the two point function. The next order contributions are coming 
from operators of dimension 8. To calculate their contribution is a quite extensive task. Also, due to 
the lack of numerical information on dimension 8 operators, one may wonder if that is worth while 
at the present stage. 
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